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Abstract. The SMAC project is a Fundamental Plane peculiar velocity 
survey of 56 clusters of galaxies to a depth of cz ~ 12000 km s-^ We 
present here some results from the analysis of the SMAC velocity field, 
I focussing on three specific features: the best-fitting bulk-flow model for 

' the SMAC data; the agreement between the observed velocity field and 

^\ [ predictions from the IRAS-PSCz redshift survey; the role of the Great 

Attractor and Shapley Concentration in generating the local flows. We 
argue that the local mass distribution, as probed by the PSCz, can fully 
account for the observed cluster velocities. 
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1. Introduction 



The 'Streaming Motions of Abell Clusters' (SMAC) project is a Fundamental 
Plane (FP) survey of ~700 early-type galaxies in 56 local rich clusters. The 
^ ' cluster sample has approximately full-sky coverage, and a limiting depth of ~ 

^ [ 12000 km s~^. Within each cluster, distances are based upon data for 4-56 E/SO 

galaxies. Data for the SMAC project is drawn from a compilation of literature 
sources and an extensive body of new observations, carefully combined into a 
homogeneous database. Further details of the sample, the observations, and 
the data compilation procedures are reported by Hudson et al. (1999) and in a 
forthcoming series of papers. 

Our (inverse) FP analysis closely follows the methods employed, for a 
smaller sample of clusters, by Hudson et al. (1997). The FP yields distances 
to a precision of ~ 20% per galaxy, so that cluster distances are determined 
with errors of 3-12%, with a median 8%. The velocity zero-point is calibrated 
by requiring that the sample exhibit no net radial inflow or outflow. The bal- 
anced sky-coverage of the sample ensures that there is little covariance between 
monopole and dipole components of the velocity field. The observed velocity 
field is presented in the upper panel of Figure 1. Lucey et al. (in this volume) 
present some example comparisons of the SMAC distances/ velocities with those 
determined from other surveys with clusters in common (eg Lauer & Postman 
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1994; Dale et al. 1999). Such comparisons generally reveal an acceptable level 
of agreement, in so much as the scatter is compatible with the quoted errors. 
However, the distance errors are such that cluster-by-cluster comparisons are 
rather crude tests for systematic effects in distance estimates. 

2. The SMAC bulk-flow 

The best-fitting bulk-flow model is the vector Vb which minimizes 

where Vi is the observed radial peculiar velocity for cluster i, whose direction 
vector is f^. The weights are assigned according to the observational peculiar 
velocity error, fTj, and a 'thermal' velocity noise cjth) here fixed at 250kms~"^. 
In what follows, we shall refer to as 'the bulk-flow', but this terminology 
should be understood as a shorthand for 'the best-fitting bulk-flow model'. As 
demonstrated by Fcldman & Watkins (1994) and further illustrated by Hudson 
et al. (this volume), the incomplete cancellation of small- and intermediate- 
scale flows causes bulk-flow measurements to depend upon the sample geometry, 
characterised by the survey window function. This effect is especially important 
when the spatial sampling is sparse, as for cluster surveys. 

For the SMAC sample, we determine the following CMB-frame bulk-flow 
solution: 

V^*"* = [-345 ± 85, +37 ± 101, -538 ± 158] km s"^ (2) 

in supergalactic cartesian co-ordinates (used throughout this paper for flow vec- 
tor components). After correcting for 'error-biasing', this vector has magnitude 
640 ± 200 kms"^ and is directed towards {l,b) = (260°, -2°). Note that the 
error ellipsoid is anisotropic, as demonstrated by Monte-Carlo simulations in 
the lower panel of Figure 1. The bulk-flow errors include a contribution from 
'system-matching' errors which introduce covariance between cluster distance 
errors (Smith et al. 1997; Hudson et al. 1997). 

The large amplitude of the observed flow is initially surprising. We have 
tested for a wide range of possible systematic effects which might affect this 
result: no single cluster, or supercluster region, is responsible for the flow; cor- 
recting FP distances for stellar population differences, based on the clusters' 
offsets from the Mg— cr relation, does not significantly affect the bulk-flow; our 
choice of Schlegel et al. (1998) extinction maps over Burstein &: Hieles (1982) 
makes little difference to the result; a 'Method IF analysis, which fits simultane- 
ously the FP parameters and flow model and is insensitive to Malmquist Bias, 
yields indistinguishable results. 

At face value, the SMAC result appears inconsistent with the bulk-flow 
solution of Lauer &i Postman (1994), whose apex lies ~ 90"^ from the SMAC flow 
direction, as well as with the non-detection of bulk-motion in the Tully-Fishcr 
survey of Dale et al. (1999, and these proceedings). In fact, the inconsistencies 
between surveys are marginal when one accounts for differences in the survey 
window functions (again sec Hudson in this volume). For the same reason, 
the apparently good agreement between the SMAC and Willick (1999) solutions 
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Figure 1. The upper panel shows the sky-projection of the SMAC 
pecuhar velocity field (CMB frame), in galactic co-ordinates. Clusters 
with positive (negative) peculiar motions are indicated by asterisks (cir- 
cles), with symbol sizes are assigned according to the magnitude of the 
cluster peculiar velocity. The S indicates the direction of the observed 
bulk- flow (ie the vector V^*) while P shows the PSCz-prediction for 

the bulk-flow (Vg^'^'^). The cloud of small crosses show the directional 
error determined from Monte-Carlo simulations. In the lower panel, we 
show the bulk-flow direction computed from 1000 Monte-Carlo simu- 
lations in which the true cluster velocities are all zero. The preference 
for directions b ~ 0°,l ~ 50° , 230° results from the relatively poor 
sampling of this axis by the cluster distribution. 



should not be over-interpreted: the surveys have very different window functions, 
and we should not expect the derived bulk-flows to agree so closely! 

3. Comparison to IRAS PSCz 

How to interpret the SMAC bulk-flow result depends upon whether the observed 
velocities of clusters (individually, or combined via the best-fitting bulk-flow 
model) can be understood as the expected response to the surrounding density 
field. Here, we compare the SMAC data to density-velocity reconstructions 
from the most extensive all-sky redshift survey of IRAS galaxies, the PSCz 
(see contributions of Saunders and Branchini in this volume). The analysis 
presented here is somewhat preliminary, and various systematic effects remain 
to be investigated. 

First, let us consider the predictions for the bulk- flow. We compute pre- 
dicted velocities for each of the SMAC clusters (assuming Pi = = 1 
for now), and substitute these predictions in place of the Vi in Equation We 
again minimize with respect to V^, using the same weights computed from 
the SMAC measurement errors and the thermal velocity noise. The best-fitting 
flow vector is what PSCz predicts for the SMAC bulk- flow, ie. given the SMAC 
errors and the SMAC sample geometry. Specifically, we fully account for any 
over-representation of clusters in (for example) the Creat Attractor direction. 
The expected bulk-flow solution for /?/ = 1 is then: 

yP^"^ = [-209, +195, -475] km s'^ . (3) 

A cursory comparison between Solutions ^ and |3| reveals that the PSCz velocity 
field for Pj = 1, given the SMAC sampling geometry and measurement errors, 
predicts a best-fitting bulk-flow similar to that observed, (including the large 
negative SGZ component). Solution ^ has magnitude 490+140 km s^^ (error- 
bias corrected) and direction {l,b) = (253°, 14°), which is only ~17° from the 
direction of the observed flow apex (see Figure 1). From Monte-Carlo simulations 
(lower panel of Figure 1), we find that such a good directional agreement arises 
by chance in < 3% of realisations. 

Alternatively, we can compare directly the predicted and observed velocity 
fields on a cluster- by-cluster basis, as in Figure 2. The fit yields a measurement 
of Pi = 0.95 + 0.19, with no individual cluster influencing the result by more 
than 10%. We can introduce greater freedom into the model velocity field by 
fitting simultaneously for Pi and a 'residual' bulk-flow vector. This residual 
bulk-flow will absorb any dipole signature not accounted for by the PSCz galaxy 
distribution (either a real streaming generated at very large depths, or a spurious 
flow resulting from systematic errors). The results of the fit are an unchanged 
value for /?/ = 0.94 + 0.25 and a residual flow vector: 

^resid ^ ^_-^gg ^ ^Q3^ ^ ^^Q^ ^ ^gQj kms"^ . (4) 

The errors are larger than previously since we are now fitting for more param- 
eters; however, the correlation of observed with predicted peculiar velocities 
remains significant at the > 3a level. That the residual bulk-flow is not signif- 
icant is a strong indication that the observed SMAC flow is indeed compatible 
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Figure 2. Cluster-by-cluster comparison of SMAC-measured peculiar 
velocities to the radial peculiar velocities predicted from IRAS-PSCz 
(for /9/=l). Symbol sizes are proportional to the weight carried in the 
fit for /?/. The lines show the best-fitting [ij and its la errors. 



with the expected response to local density fluctuations; and suggests that the 
observed flow does not result primarily from systematic distance errors. The 
favoured value for f3i is higher than many estimates, but is consistent with other 
determinations based on PSCz, which yield /?/ = 0.6 — 0.7 (see Branchini, this 
volume). Note that there is some covariance between /?/ and V^''*'^, such that 
flxing (3i = 0.7 results in a solution with larger residual flow, albeit significant 
at only the l.So" level. 

A final caveat should be noted: the preliminary analysis presented here 
is based on a Method I approach and is therefore affected by Malmquist Bias. 
Homogeneous Malmquist effects have been corrected for and the effects of Inho- 
mogeneous Malmquist Bias (1MB) are suppressed by the use of a cluster sample 
with relatively small distance errors per 'object'. Nonetheless, it may be ex- 
pected that the 'spurious infalP patterns produced by 1MB act to artificially 
infiate /?/. A crude attempt to judge the effect can be made by removing the 
12 clusters with |vp5Cz| > 500kms~^, which lie mostly in infall regions where 
1MB will be most severe. The result has much greater uncertainties of course, 
since we have removed the clusters which contribute the greatest signal-to-noise: 
we obtain /?/ = 0.6 it 0.4, with no measurable residual bulk-flow (amplitude 
100 lb 200kms~^). In a future analysis of the SMAC-PSCz comparison, we will 
employ a Method II approach which is free from Malmquist bias effects. 
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4. The Great Attractor and Shapley Concentration 

The SMAC sample includes 10 clusters within 15° of the Shapley / Great At- 
tractor (GA) direction (visible in Figure 1 as the concentration of points at 
I = 315°, 6 = 30°). Since this region is of some historical interest, we compare in 
Figure 3 the observed velocities with some (over-) simplistic toy-models for the 
dynamics of the region. The models are based upon the simple spherical attrac- 
tors of Faber &: Burstein (1988). We consider the effect of two such structures, 
one at the position of the 'traditional' GA (at a distance of 4500 kms~^) and 
a second centered on the Shapley core region (at 14500 km s^^). The models 
are normalised to generate the GA-directed component of the Local Group's 
velocity in the CMB frame, viz. 500 km s"-*^. 

The leftmost panel of Figure 3 demonstrates that a pure GA infall model is 
a poor fit to the SMAC data, especially in the immediate GA background, where 
no 'far-side infall' is observed. Equally, the extreme Shapley infall model dramat- 
ically over-predicts velocities beyond ~10000kms^^ (middle panel). Allowing 
for contributions from both attractors, and fitting for their relative amplitudes 
as the single free parameter, we can obtain an acceptable fit to the data, as 
shown in the rightmost panel. The best-fitting model has Shapley and the GA 
each generating 50 it 10% of the Local Group velocity in this direction, pro- 
portions which accord with Hudson's (1994) conclusions based on the Mark II 
dataset. 

Here, as in the case of 'bulk-flow' wc caution against drawing quantitative 
conclusions from analyses based on unrealistically simple models of the velocity 
field. However, the above fits draw attention to the qualitative behaviour of 
the SMAC velocity data in the GA / Shapley direction. The absence of far-side 
infall in the SMAC data (here visualised as a retardation of one infall pattern by 
another more distant infall structure) is also apparent in the PSCz maps, which 
reveal a wealth of structure along the GA-Shapley axis. 

5. Conclusions 

From our analyses of the SMAC velocity field, as presented in the preceeding 
sections, we draw the following conclusions: 

1. The best-fitting bulk-flow model has amplitude 640 ± 200 kms"^, towards 

{l,b) = (260°, -2°). 

2. The direction of the observed bulk-flow is within '-^15° of the direction 
predicted from the IRAS-PSCz redshift survey for the SMAC sample, 
while the amplitude can be matched if /?/ ^ 1. 

3. Comparing observed and predicted velocities cluster-by-cluster yields /3 = 
0.95 ± 0.25, compatible with other analyses based on the PSCz. 

4. The observed velocity field is fully accounted for by the local density field, 
with no recourse to residual flow generated beyond the limits of PSCz, or 
to 'spurious' flows associated with systematic errors. 
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Figure 3. SMAC peculiar velocities in the Great Attractor / Shapley 
direction. The data are identical in each panel. The models shown are: 
a spherically symmetric Great Attractor model (left); a similar model 
attractor at Shapley (middle); the best fitting 'two-attractor' model, 
with errors (right). All models are normalized to generate 500kms~^ 
at the Local Group. 

5. The Shapley Concentration may be responsible for a significant fraction of 
the infall traditionally associated with the Great Attractor. 
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